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NOTATION 


dr 

d1 

d 

E 

v 

R 

h 

L 

Z 

e 


V  V  Yx0 


displacement  from  rest  (generic  form) 
initial  imperfection  (generic  form) 
axisymmetric  displacement  (generic  form) 

Young's  Modulus 
Poisson's  ratio 
shell  radius 
shell  thickness 
shell  section  length 

distance  from  shell  mid-surface  (positive  inwards) 

distance  from  shell  mid-surface  to  frame  centroid 
(positive  for  internal  frames) 

circumferential,  axial  and  shear  strains  respectively 


x 


e 


i*  i’ 


frame  area 

axial  position;  as  a  subscript  denotes  differentiation 

with  respect  to  x  when  used  in  conjunction  with  u,  v,  w, 

i 

w 

circumferential  co-ordinate;  as  a  subscript  denotes 
differentiation  with  respect  to  0  when  used  in  con¬ 
junction  with  u,  v,  w,  w1 

eigen  vector  amplitudes  associated  with  I1*1  basis 
function  in  the  expansion  of  axial,  circumferential  and 
radial  bifurcation  displacements,  respectively 


u(x,6),  v(x,e),  w(x,9)  axial,  circumferential  and  radial  bifurcation  displace¬ 
ments,  respectively,  (x,  and  6  often  omitted), 
w  positive  inwards 

wX(x,9)  the  radial  out-of-circularity  at  x  and  0.  (x  and  0 

often  omitted) 


i 

w 

X 


derivatives  of  w 


with  respect  to  x  and  0,  respectively 


u(x) 
W(x) 
w  (x) 

X 


axial  pre-buckling  displacement 
radial  pre-buckling  displacement 
pre-buckling  rotation 
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upper  limit  to  the  number  of  basis  functions  required  in 
the  expansion  of  bifurcation  displacements 

number  of  circumferential  waves 

maximum  departure  from  mean  circle 


displacement  vector 

lateral  pressure  (positive  external) 

axial  pressure  (positive  external) 

wave  number  of  imperfection  shape,  eoc  sin  I^itx/L 


number  of  frames 


number  of  frame  spaces 


frame  spacing 


width  of  faying  flange ,  or  web  thickness 
centroidal  radius  of  frame 


second  moment  of  area  of  frame ,  for  in-plane  bending 


frame  area 


in  the  Kendrick  Frame 


,  .  ,  .  expressions  (Appendix) 

second  moment  of  area  )  * 
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CLOSED  FORM  SOLUTION  FOR  THE  BENDING  STRESSES  IN 
RING-STIFFENED  SHELLS  WITH  SHAPE  IMPERFECTIONS 


INTRODUCTION 


An  important  consideration  in  the  design  of  externally  pressurised 
cylindrical  shells  is  the  allowance  for  the  shape  imperfections  which  inevitably 
occur  during  manufacture.  Components  of  out-of-circularity  which  approximate  to 
critical  buckling  mode  shapes  are  considerably  exaggerated  as  the  applied 
pressure  approaches  the  corresponding  buckling  pressure.  Frame  stresses  may 
exceed  the  material  yield  strength,  precipitating  failure  by  one  of  several 
mechanisms. 

2.  Many  methods  are  available  for  finding  the  buckling  pressures  of  the  perfect 
ring-stiffened  shell,  and  while  frame  stresses  may  be  found  with  simplifying 
assumptions  concerning  the  way  mode  shapes  will  develop  (1,  2,  3)*,  comparatively 
few  analyses  exist  which  allow  a  completely  satisfactory  direct  solution  to  the 
bending  stresses. 

3.  The  purpose  of  the  analysis  presented  in  this  paper  is  to  provide  a  closed 
form  solution  for  the  displacements  and  frame  bending  stresses  in  a  nearly 
circular,  uniformly  framed,  ring-stiffened  cylindrical  shell.  The  method  is 
based  rigorously  on  thin-shell  theory,  but  is  nevertheless  simple  enough  to 
allow  calculations  to  be  made  on  a  programmable  pocket  calculator. 

THEORY 


4.  Kendrick  has  shown  that  the  strains  in  a  slightly  non-circular  shell  may  be 
written  as  follows  (1): 


=  (vQ  -  w)/R  + 


(w. 


+  u„ 


2wve)/2R* 


(w0\j0  -  w*w)/R2 


e 


x 


u  t  (l/2)(w  2 

X  X 


2 

+  v 


X 


)  + 


y  -  u„/R  +  v 
x  e  x 


M 


The  strains  contain  quadratic  products  of  the  displacements,  u,  v,  w  and  the 

imperfections  w1.  From  equations  such  as  [Y]  ,  a  set  of  linear  equations  may  be 
derived,  the  solution  of  which  yields  the  displacements  and  stresses  in  the 
shell.  An  outline  of  the  derivation  is  given  below  in  general  terms,  applying 
to  any  set  of  strain  expressions.  A  specific  derivation,  using  the  original 
Kendrick  expressions  (1)  is  given  in  the  Appendix. 


5.  An  expression  for  the  strain  of  shell  and  frames  is  written  in  terms  of  the 

/  r  i  v 

displacements  from  rest  (d  )  and  the  initial  imperfection  (d  ).  Parameters  d 

and  da  represent  in  generic  form  terms  of  the  following  type : 


*(  )  =  References  on  page  16 
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r  r  r  r 

d  :  w  ,  v  ,  u  and  their  derivatives  with  respect  to  axial  position  x 
and  circumferential  position  6 

,i  i  i  i 

d  :  w  ,  w  0,  w  x 

where  w,  v  and  u  are  respectively  radial,  circumferential  and  axial  displace¬ 
ments. 

2 

•  r*  r*  r  i 

The  strain  expression  contains  terms  proportional  to  d  ,  (d  )  and  d  d  ,  the 

exact  form  of  the  expression  depending  to  some  extent  on  the  basic  formulation 

used:  the  Novazhilov  formulation  uses  a  Cartesian  co-ordinate  system;  the 

method  used  by  Kendrick  uses  cylindrical  co-ordinates. 

6.  The  total  potential  energy  is  given  by  the  following  expression 


U 


T 


U 

s 


+  U, 


W 


where  the  contributions  from  shell  (Ug),  frames  (U^.)  and  external  load  (W)  may 
be  expressed  as 


2ir 


u  (x,  0)  dxd0 
s 


2ir 


I  u  (x,  0)  d6 
frames 


2ir 


t  L 


w  (x,  0)  dxd0 


M 


The  shell  and  frame  energy  densities  (u  and  u^  respectively)  contain  the  follow- 

•  •  r 

ing  types  of  terms  up  to  cubic  order  in  d  : 


(dr) 


(dr) 


<dr>V 


Cdr)V 


The  load  potential  term  has  terms  in: 
dr  (dr)  and  dr  d1 

7.  Total  displacements  (d  )  are  expressed  as  the  sum  of  primary  displacements 
(d)  and  bifurcation  or  buckling  displacements  (d): 

dr  =  d  t  d  [5] 

The  following  terms,  of  less  than  cubic  order  in  d  or  d,  are  then  present  in  U^: 
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d  d 


d2 

dd 

d2 

d3 

d2d 

dd2 

d3 

d'd2 

dXdd 

d\l2 

d^d3 

d^2d 

dXdd2 

d\i3 

1 

1 

4 


8.  The  primary  displacements  (d)  at  any  desired  pressure  p  are  evaluated 
first,  including  if  necessary  the  effects  due  to  an  axisymmetric  imperfection. 

It  is  assumed  that  non  axisymmetric  imperfections  are  sufficiently  small  so  that 
d  may  be  considered  to  be  axisymmetric.  The  bifurcation  displacements  d  are 
found  by  minimising  U^.  Bifurcation  displacements  are  assumed  to  be  sufficiently 

small  to  enable  coupling  between  different  lobar  modes  to  be  ignored.  The 
bifurcation  mode  shape  is  considered  therefore  to  be  a  pure  Fourier  mode.  With 
the  simplifying  assumptions  concerning  pre-buckling  and  buckling  displacements, 
only  five  sets  of  terms  need  be  considered: 

From  the  strain  energy  densities: 


1.  d2  (Eh/R) 

2.  dd  (Eh/R2) 

3.  dd1d  (Eh/R2) 

4.  d\i2  (Eh/R2) 

and  from  the  external  work: 

2 

5.  pd 

where  p  is  the  applied  pressure. 


The  typical  order  of  magnitude  of  coefficients  corresponding  to  each  type  is 
indicated  in  brackets.  Terms  1  and  4  are  independent  of  external  load  (p). 
Terms  2,  3  and  5  are  directly  dependent  on  p.  The  coefficients  of  type  4  terms 
are  such  that  the  contributions  of  the  terms  are  smaller  than  the  contributions 

from  type  1  terms,  by  the  ratio  d1/R.  The  effect  of  type  4  terms  is  to  couple 
different  lobar  modes  for  large  displacements.  Provided  the  level  of  out-of- 
circularity  is  small  in  relation  to  the  radius,  and  the  gradient  of  the 
imperfection  surface  not  too  rapidly  varying  in  the  x-direction,  type  4  terms 
may  be  ignored,  maintaining  compatibility  with  the  assumption  that  the  develop¬ 
ment  of  different  lobar  modes  occurs  independently. 

9.  Bifurcation  displacements  are  represented  in  terms  of  Fourier  functions. 
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u  lx,  0)  =  cos  n0  I  u  cos  Iirx/L 

1=1  L 

NI 

v  (x,  6)  =  sin  n0  I  v  sin  Iirx/L  [6] 

1=1 

NI 

w  (x,  0)  =  cos  n0  I  w  sin  Iirx/L 

1=1  1 

On  minimisation  of  with  respect  to  the  amplitudes  u^. ,  v^ ,  w^  ,  for  I  =  1  to  N 
and  scaling  to  a  non-dimensional  form ,  the  following  matrix  problem  arises 

M  x  =  £  [7] 

where  M  =  A  +  4/B 

4/  =  pR  (1  -  v2)  /  Eh  [8] 


Matrix  A  contains  coefficients  of  type  1  terms  (d  ).  Matrix  B  comprises 
coefficients  of  terms  from  types  2  and  5.  The  vector  c_  derives  from  group  3 
terms,  and  _x  is  the  displacement  vector  whose  terms  are  the  amplitudes  of  the 
Fourier  functions  at  any  pressure  p: 


T 

x 


(U1  V1  W1 


*1  VI  WI 


10.  For  the  perfectly  axisymmetric  structure  (£  =  0)  the  matrix  problem  of 
equation  Jj7]  reverts  to  the  standard  buckling  eigen  value  problem.  Matrices  A 
and  B  are  identical  to  the  buckling  matrices. 

11.  Prebuckling  distributions  required  for  the  B  matrix  and  c_  vector  have,  in 
this  paper,  been  calculated  as  follows.  The  axisymmetric  radial  deformation  of 
a  frame  in  a  uniformly  framed,  ring-stiffened  cylinder,  under  unit  hydrostatic 
load  is  given  precisely  by  the  following  expression  (4,  5) 

w  =  R2  (1  -  v/2)  /  {  Eh  [l  +  A  /  (bh  +  2Nh/ct)  ]  }  [lo] 

where 

(R/R  )2 
g 

'  -  v2>?  /  /Rh 

N  =  (cosh  aL*  -  cos  aL' )  /  (sinh  oL'  +  sin  aL’) 

L'  =  L  -  b 
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To  conserve  radial  load,  the  net  circumferential  stress  in  the  shell  must 
satisfy  the  following  equation: 

W  =  “•f  -  /  Rg  & 

'*  _  sh 

The  longitudinal  pre-stress  ox  is  given  by 

ixSh  =  -R/2h  [l 

TJie  axisymmetric  circumferential  strain  (e  s^)  and  the  axial  strain  (u  S^)  in 
the  shell  are:  *  K 


-  sh  ■ 

ee  =  E 


-1  [1eSh  -  vR/2hj 
-1  [v  SeSh  -  R/2h] 


uxSh  =  E_1  [v  a0Sh  -  R/2h  [l4] 

for  the  chosen  sign  convention,  pressure  positive  external,  displacements 
positive  inwards. 

12.  The  matrix  M  may  now  be  defined  very  simply  (details  in  Appendix).  Shell 
and  frame  components  are  considered  separately, 

M  =  M  +  M,.  . 
s  f 

For  a  ring- stiffened  cylinder  of  uniform  shell  thickness  the  shell  matrices  are 
banded  tri-diagonally , 

M  (IU  ,  IU)  =  m  UU  +  I2m  UU 
s  u  z 

l iv 

M  (IU  ,  IV)  =  Im. 

s  1 

It,,  t,,\  t  UW  ,3  UW 

M  (IU  ,  IW)  =  Im  +  I  m_ 
s  1  3 

D-d 

,  vv  2  w 

M  (IV  ,  IV)  =  m  +  I  m 

s  o  2 

M  (IV  ,  IW)  =  m  W  +  I2n  ™ 

s  o  2 

w  /T,.  TT.7\  ww  _2  ww  _4  ww 

M  (IW  ,  IW)  =  m  +  I  m„  +  I  m. 

S  0  2  4 


where 


=  31-2 


=  31-1 


I  =  1,  2, 


UNLIMITED 
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to  be  compatible  with  the  displacement  vector  definition,  equation  [9]. 
uu 

Expressions  for  m^  etc  specific  to  the  Kendrick  formu1  - ti«_>n  are  given  in  the 

Appendix. 

13.  Satisfactory  solution  of  overall  buckling  problems  may  be  obtained  if  out- 
of-plane  bending  stiffness  of  frames  is  ignored  in  both  A  and  B  matrices  (3). 
The  frame  terms  of  the  M  matrix  then  contain  neither  x-derivat ives  nor  terms  in 
axial  displacement  u. 


M  (IV  ,  JV)  =  fV  S(I,J) 


M  (IV  ,  JW) 


VW 

f  S(I,J) 


M  (IW  ,  JW)  =  f  S(I,J) 


where 


S(I,J)  =  (2/N  )  Z  sin  (Iirx./L)  sin  (Jttx./L) 

S  i=l 


[17] 


[18] 


14.  The  c_  vector  may  be  similarly  decomposed  into  shell  and  frame  terms 


c  =  c  +  C-j. 

—  —s  — f 


If  the  imperfection  is  of  the  form  e_^  cos  n9  sin  (I^irx/L)  the  shell  components 

of  the  c_  vector  may  be  written  as  follows  (equations  [1 3]  and  [23]  of  the 
Appendix): 


(IU)  = 

T  UX 

1  C16II 

(IV)  = 

cV5 
o  II 

oc 

(IW)  = 

[  W  .  ' 

c  +  c 

oc 


[19] 


II 


oc 


The  frame  terms  are  given  by 


c£  (IV)  =  cl  c(I) 
f  f 


c  (IW)  =  c*  c(I) 


[20] 


where 


N’f 

c(I)  =  (2/N  )  Z  sin  (Iirx./L)  sin  (I  ttx./L) 


i  =  l 


oc  1 


[21] 
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15.  For  a  uniformly  framed  structure  (identical,  uniformly  spaced  frame* 
sinusoidal  components  obeying  the  following  relationship  for  some  integral  value 
of  K  are  coupled  (3): 

I  =  (^)  {jl  -  (  -l)1^  [<I>  -DNS  +  kJ  ±  |l  +  (  -l)1^  [I'Ng  -  kJ  } 


for  I*  =  1,  2,  3, 


Functions  satisfying  equation  [22]  with  either  "  "  or  in  the  centre  provide 

a  complete  orthogonal  set.  Selecting  the  option,  it  may  be  shown  straight¬ 
forwardly  that 

S(I,J)  =  1  [23] 

for  all  I,J  satisfying  equation  [22].  Only  the  long  wavelength  shapes  will  be 
required , 


eQc  sin  (Kirx/L)  for  K  =  1,  2,  3, 


representing  the  fundamental  modal  components  of  equation  [22].  If  I  =  K, 

K  <  N  ,  then  oc 

c(D  =  1  [24] 

for  all  I  satisfying  equation  [22].  It  follows  immediately  that  the  fundamental 
matrix,  equation  [7] ,  may  be  stated  in  the  following  form: 

M  (IU,IV)u  +  M  (IV,IV)v  t  M  (IV,IW)wT  +  f  VEvT  +  fVWZwT 
S  IS  X  S  X  X  X 

2  co5iioc  *  cf 

M  (IU,IW)u  +  M  (IV,IW)v  +  M  ( IW,IW)wT  +  fVWIvT  +  fWEwT 
S  X  s  xs  X  X  X 

=  (c^c«I2)6ii  +  c”  [26] 

oc 

M  (IU,IU)uT  +  M  ( IU ,1V) v  +  M  (IU,IW)wT  =  I  c*6Tr  [27] 

s  is  I  s  I  III  L--J 

oc 

the  summations  being  performed  for  a  set  of  I  given  by  [22]  . 

16.  From  equation  [27]  ,  can  be  written  explicitly  in  terms  of  v^  and  w^ , 
allowing  the  matrix  dimension  to  be  reduced  by  a  factor  of  2/3: 

S  (Dvj  +  S  (I)Wj  =  c  (I)  -  f  EVj  -  f  Ewj 


SVW(I)vI  t  SW(I)wI  =  cW(I)  -  fVWEvI  -  fWEWj 
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where 


sv(i) 


M  (IV, IV)  -  M  (IU.IV)  /  M  (IU,IU) 
s  s  s 


SWII> 


M  (IV, IW)  -  M  (IU.IV)  M  (IU,IW)  /  M  (IU,IU) 
s  s  s  s 


SW(I) 


Ms(IW,IW)  -  M^dU.IW)  /  Ms(IU,IU) 


cV(I) 


V  V 

C,  +  c  , 

f  sh 


c\l) 


w  w 

C.P  +  c  . 
f  sh 


c  h  -  5tt 
sh  II 

oc 

V  _  u 

[fo  -  1  ci 

Ms(IU,IV) 

/  MS(IU,IU) 

W  _  r 

Csh  “  6II 

oc 

w  w  t2 

[fo  +  °2  1 

-  I  c'f  M 

1  s 

(IU,IW)  /  Ms 

17.  Equations  [28] 

show  that 

Vj  and  w^. , 

for  any  I , 

of  summations  over  1 

and  w  : 

v( I )  =  sW(I) 

[cv(D  -  5 

]  -  •v"«) 

Ewu)  -  3 

w(I)  =  sv(I) 

Q“(I)  -  £ 

]  -  s™U) 

[«*»>  -  3 

where 

-V  v  .vw  _w 

n  =  f  E  +  f  e 

C  =  fW  IV  +  fW  EW 


v  w 

I  and  E  denoting  Ev^  and  Ew^  respectively. 


sV(I)  =  SV(I)  /  D( I ) 


sW(I)  =  S™(I)  /  D( I ) 


sW(I)  =  SW(I)  /  D ( I ) 


D( I )  =  SV(I)  SW(I)  -  S™(I)2 
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18.  Substituting  equations  [32]  into  equations  [_28]  ,  it  may  be  seen  that 

rv  «  {fl  ♦  5  fw-i  fw]  c.  -  [s  fw-i  fw]cJ/ 

lL  v  vw  J  1  [_w  vw  J  2  / 

{ [l  +  i  fV  -  i  f™]  fl  +  l  fW  -  i  fW|  -  [35] 

1  j_  w  vwj[_v  VW  J 

[1  fw-i  fv]fi  fw-i  /h 

[_  v  vwj[_w  VW  _|  J 

zw  =  [c  -  (1  +  s  fv  -  i  fW)  1^1  /  (S  f w  -  i  fw) 
l_l  w  vw  J  /  w  vw 

where 

s  =  Isv(I)  s  =  Isw(I)  i  =  Esw(I) 

V  vw  w 


[38] 

19.  Equations  [32]  to  [38]  represent  the  closed  form  solution  for  the  circum¬ 
ferential  (v)  and  radial  (w)  displacements  of  a  hydrostatically  loaded,  slightly 
non-circular,  uniformly  framed,  ring-stiffened  cylinder,  with  shape  imperfection 
eQc  sin  (Kttx/L)  cos  (n0).  Equations  [35]  to  [38]  may  be  evaluated  directly  in 

terms  of  known  shell  and  frame  coefficients  (equations  [21]  and  [25]  of 
Appendix). 

20.  Equations  [32]  give  v(I)  and  w(I)  from  which  the  shell  displacements  at 
any  position  x  using  equations  [6]  and  [22]  .  In  this  derivation  the  central 
sign  of  equation  [22]  has  been  taken  to  be  minus.  If  the  ’plus'  option  is  used 
the  final  expressions  take  on  the  slightly  modified,  but  numerically  equivalent, 
form  given  in  (3). 

21.  The  stresses,  at  position  z  on  a  frame  section,  are  given  by  the  following 
expression 

cr(z)  =  a  (z)  +  a,  (z)  [39] 

ax  b  u  J 

where  the  axisymmetric  stress  a  (z)  is  given  by 

clX 

a  (z)  =  E  wc  /  (R  -  S.z)  ,  [40] 

ax  f  /  i  u  J 

and  the  bending  stress  a^Cz)  Siven  by 

o^Cz)  =  (E/R)  (n  v(x)  -  w(x)  +  (z/R)  (n2  -  1)  w(xf]  [4l] 
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V  -  W  ^  W,,  V  vw,,,  w 

s  c  -  s  c  +  s  (1)  c  -  s  (1)  c 
w  f  vw  f  sh  sh 


w  -  V  v.,Nw  VW, , ,  V 

s  c  -  s  c  +  s  (1)  c  -  s  (1)  c 
v  f  vw  f  sh  sh 
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NUMERICAL  EXAMPLES 


22.  Displacements  and  stresses  have  been  evaluated  on  the  two  rather  different 
structures  given  in  Table  1. 

23.  Displacements  and  stresses  in  the  closely  frame-spaced  submersible 
structure  are  given  in  Tables  2-4.  The  predictions  of  closed  form  and  matrix 
method  are  almost  identical.  The  closed  form  method  is  particularly  useful  for 
convergence  tests  (Table  4). 

24.  In  the  widely  frame-spaced  VAC  type  structure  severe  mixed-mode  effect 
occurs  for  n  >  2.  The  closed  form  method  provides  therefore  a  very  simple 
method  of  allowing  for  mixed  mode  effect  in  ring-stiffened  cylinders.  Displace 
ments  and  stresses  predicted  by  closed  form  and  matrix  methods  agree  almost 
exactly  (Figure  1,  Table  5). 

CONCLUSIONS 


25.  By  means  of  analytical  simplification,  it  is  possible  to  derive  rigorously 
from  thin-shell-theory ,  a  closed  form  expression  for  the  bending  stresses  in  a 
uniformly  framed  ring-stiffened  cylinder  with  non-axisymmetric  shape  imperfec¬ 
tion. 


26.  The  analytic  method  has  been  implemented  on  a  programmable  pocket  calcu¬ 
lator  and  results  are  shown  to  agree  with  the  predictions  of  more  usual  matrix 
methods  which  require  quite  large  computers  for  accurate  solution. 

27.  It  is  felt  that  the  'closed  form'  method  provides  a  useful  means  of  check¬ 
ing  the  prediction  of: 

a.  more  sophisticated  methods,  such  as  finite  difference  and  finite 
element ,  and 

b.  approximate  methods  used  in  design  codes,  such  as  the  effective 
breadth  approach  of  BS5500. 

The  method  is  particularly  appropriate  for  carrying  out  convergence  studies  to 
assess  the  fine-ness  of  mesh  required  in  more  sophisticated  analyses. 
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tabu:  i 


CASE  1  -  PV4  SUBMERSIBU 
19  equispaced  frames 


R  =  2000  mm 

h  =  17 . 879  mm 

L  -  12000  mm 

E  =  207000  N/mm2 

v  =  0.3 


d  =  180.92  mm 

h  =  9.046  mm 

w 

f  :  90.46  mm 

-  15.07  nan 

INTERNAL 


CASE  2  -  PV4  VAC  TYPE 
1  central  frame 

P  =  2500  mm 

fc  :  20  mm 

L  =  14000  mm 

E  z  207000  N/mm2 
v  =  0.3 


d  =  192  mr. 

h  =  12  mm 

w 

i  -  0 

hf  =  o 

INTERNAL 


R  shell  radius 

h  shell  thickness 

L  shell  length 

E  Young's  modulus 

v  Poisson's  ratio 


d  web  depth 
h^  web  thickness 
f  flange  width 

flange  thickness 
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PV4  Submersible  n  =  2  displacements  and  bending  stresses  at  3  N/mrn  for  shape  20  cos  20  sin  2irx/L. 
Comparison  of  results  from  the  closed  form  method  (M  =  5)  and  a  matrix  method  ( 30  x  30  matrix, 

IFAC  =  0). 
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TABLE  4 


M 

N 

V 

W 

«>fl 

c  fl 

TOT 

1 

6 

3.854 

8,060 

220 

417.0 

2 

12 

4.080 

8.653 

224 

421.5 

3 

•  ' 

4.106 

8.723 

225 

422.0 

5 

" 

30 

4.116 

8.749 

225 

423.0 

PV4  Submersible  n  =  2,  shape  20  cos  20  sin  nx/L. 

Convergence  test:  M  convergence  parameter, 

N  equivalent  (IFAC  =0)  matrix  size. 

2 

Pressure  3  N/mm  . 
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TABLE  5 


FRAME  DISPLACEMENTS  -  PV4  VAC  TYPE  STRUCTURE 


pressure  =  0.3  N/mm 


e  =  25  mm 
oc 


I  =1 
oc 


n  =  2 


Circumferential  Displacements  (v) 


Closed  Form 
Matrix 


1.3725 

1.3724 


Radial  Displacements  (w) 


Closed  Form 
Matrix 


2.7460 

2.7449 


5.0266 

5.0292 

15.470 

15.475 


2.4920 

2.5007 

10.412 

10.444 


1.1397 

1.1502 

6.0950 

6.1474 


0.44280 

0.44616 

2.9730 

2.9940 
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APPENDIX 


In  equations  [l]  e  ,  e  ,  y  .  are  respectively  the  circumferential,  axial 

U  X  Xo 

and  shear  strains.  The  displacements  w,  v,  etc  are  measured  from  rest,  the 
superscript  "r"  being  omitted  for  convenience.  The  extensional  shell  energy 
is  then  given  by  the  standard  plane-stress,  thin  shell  expression: 


U 

e 


2ir 


[ERh/2  (1  -  v2)J  { 


%  +  ex2  +  2Ve6ex  + 


(1/2) (1  -  v)y  a2}  d6dx 

Xt7 


CAl] 


2.  The  following  expression  for  the  shell  bending  strain  energy  was  used  by 
Kendrick  ( 6 ) : 


2n  ,L 


ERhk/2  (1 


-  v2)]  {R2w 


2  2  2 
+  (w„_  +  w)  /R  + 
xx  00 


o  '  o 


(1/2)<1  -  v)(w  -  u./a)2  +  (3/2K1  -  v)(v  +  w  Q)2  +  Fa2] 
X0  0  X  X0  *-  J 


2vw  (w0.  +  v.)  +  2Ru  w  }  dxd0 
xx  00  0  x  xxJ 


2  2 

where  k  =  h  /12R  .  The  shell  energy  contribution  Ug  =  Ue  +  Ufc.  Frame  exten¬ 
sional  strains  are  expressed  in  terms  of  the  mid-shell  strain  e  and  the  change 

U 

of  curvature  due  to  in-plane  bending,  giving  the  following  expression  for  the 
extensional  frame  energy  F^, 


•2ir 

{eg  ~  (wee  +  w)(e/R2)}2  d0  [A3] 

o 

3.  Frame  bending  energy  is  given  by: 

■2tt 

(w00  +  w)2  d0  [a 4] 

o 

4.  The  total  frame  energy  =  Fg  +  F^.  The  work  done  by  external  load  on  an 
initially  imperfect  cylinder  is  found  to  be  given  by  the  following  expression: 


Z  (EI/2R  ) 
r=l 


f 

F  =  I  (EAR/2) 
r=l 
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FT 


i 


t 2tt 


W  =  (p/2) 


UNLIMITED 


|2Rw  -  2w1w  -  w^  +  2R  (w1  +  w)(u  +  v  /a) 

x  9 


Ru^Vg}  d6dx  +  (PR  /2) 


2ti 


-  u  dxde 
x 


5.  In  all  the  expressions  given  above  the  superscript  ”r"  denoting  displace¬ 
ments  from  rest  has  been  omitted.  On  making  the  transformation: 

ur(x,0)  =  u(x)  +  u(x ,9 ) 

vr(x,0)  =  v(x ,0 )  [A6] 

wr(x,0)  =  w(x)  +  w(x,6) 

the  total  potential  energy  U„  is  found  to  be  given  by  the  following  expression: 


L’a  +  UB  +  Uc 


where 


(27T 

•L 

j^EhR/2  (1  -  v2) 

o  J 

o 

{u  +  (vq  _  w)  /R  +  2vu  (v  -  w)/R  + 
x  0  x  0 


(1/2)(1  -  v)(u  /R  +  v  T  +  k 
0  x 


*V +  (w09 +  w)2/r2  + 


(1/2K1  -  v)(w  _  -  u  /R)2  +  ( 3/2)  (1  -  v ) ( v  +  w  )2  t 
X0  0  X  X0 


2vw  (w„„  +  v  )  +  2Ru  w 
xx  90  9  x  xx 


Nf 

}  d6dx  +  E  (EAR/2) 
r=l 


[A8] 


2tt 


(w  +  w ) (e/R  )  +  (w  -  v  ) /R|  + 

0  V  0 


(I/R4A)(wee  +  w)2}  d0 


-  23  - 


UNLIMITED 


2 IT 


2  2  2 

,e  -  ue  -  w0  t  4„v0  - 


J^Eh/2R2(l  -  v2^j  {w  ^2v0 

)  -  2Rvu  vj  +  u  |R3(v2+w2) 

X  _ej  X  LI  X  X 

w  2R3u  w  -  2vR2w  (w  -  v  )  \  d0dx  +  [A91 

x  J_  x  x  x  6_J  ‘  L  J 


d2  .  2  2 
R  v  (w  +  v 

X  X 


2  2 

Rv  ( u  +  w.  -  2wv  ) |  t 
yy  u 


Nf 

I  (EA/2R2 )w 
r=l 


2ir 


{-(1  +  e/R)u0  -  2v0  -  (1  +  e/R)w02  + 


4(1  +  e/R)wv  +  (2e/R)w  v  V  de  +  (p/2) 
0  00  0 


jw  -  2w  (vA  +  Ru  )  +  Ru  v  l  d0dx 
1  6  x  x  0J 


2tt 


o  *  o 


2tt 


o  4 


^Eh/2R2(l  -  v27j  fw  J^vRu^w*  -  2v^w^  +  4WW*  -  2w0w0^ 

2R2vw  w  1  +  u  -2RVWW1  +  2Rvw.vj  ^  +  2R3w  w  *  + 

x  x_J  x  |_  0  0  x  xj 

w  2R2v(v„  -  w)w  1  +  2R3u  w  1  \  d©dx  + 
x  l_  8  *  x  x  J 


2  - 

1  (EA/2R  )w 
r=l 


2tt 


|-2v  w1  +  4(1  +  e/R)ww1  - 


m 


2(1  +  e/R)w  w  1  t  (2e/R)  vrw1}  d0  + 
0  0  0  0 


2tt 


w  (w  -  v„  -  Ru  )  d0dx 
0  x 


Terms  proportional  to  the  end  loads  (P)  do  not  appear  explicitly  in  the  final 

expressions.  Effects  due  to  end  loads  are  transferred  through  w,  u  and  w  . 

x  x 
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6.  On  minimization  of  the  energy  expressions,  the  A  and  B  matrices  and  the 
c-vector  may  be  written  as  follows: 

9  2  2 

A  (IU,IU)  =  n  (1  +  k)(l  -  v)/2  +  I  X 
s 

A  (IU.IV)  =  -nIX(l  +  v)/2 
s 

2  3 

A  (IU,IW)  =  vIX  -n  (1  -  v)kIX/2  +  k(lX) 

s  M 

A  (IV, IV)  =  n2  +  (1  -  v)(l  +  3k)(IX)2/2 
s 

A  (IV, IW)  =  -nfl  +  (3  -  v)k  (IX)2/2j 

A  (IW,IW)  =  (It  k)  +  n2k  [n2  -  2  +  2 ( I A >^]  +  k(IX) 
s 


B  (IU,IU) 
s 

=  -n2b 

s  S0 

B  (IU,IV) 
s 

—  ch 

=  -nIX/2  +  n  b  vIXe 

S  0 

B  (IU.IW) 

=  IX 

S 

B  (IV, IV) 
s 

=  b  s 
s  x 

,  T,  ,2  .  2  -  sh 

(IX)  -  2n  b  e 
s  0 

jAi:j 

B  (IV, IW) 
s 

=  -nil 

,,  -  sh  ,  -  sh] 

-  2fc  e„  t  b  vu 

s  6  s  x  J 

B  (IW.IW) 
s 

=  1th 

s  -  n2  sj 
s  X  0J 

C  (IU)  = 
s 

tie  IX 
oc 

1  t  b  ve  S^j  6 

L  =  6  -1  “oc 

C  (IV)  = 
s 

' 

-tie  n 
oc 

E  *  v“] 

1^13} 

C  (IW)  = 
s 

*eoc  E 

t  b  {(2  -  n2  )e  sh  +  v(n":  -  l)u  sh  t  s  (IX)2} 

S  0  x  *  _ 

6II 

J  oc 
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A  (IV,JV)  =  n  f  S(I,J) 


Af(IV,JW) 

=  A  (IW.JV) 

=  -nf JT  -  (n2  -  l)e^j  S(I,J) 

Af(IW,JW) 

=  {f|7l  +  e 

2  2  2  2  ' 
)  -  2n  e  +  n  (n  -  2)e‘ 

r  r  j 

[]  +  S(n 

B  CIV.JV) 

=  2n2bf  S(I 

,J) 

Bf ( IV , JW) 

=  B  (IW.JV) 

=  -nb  1 2 ( 1  +  e  )  -  n2e 
f  [_  r  r_ 

S(I,J) 

B  ( IW,JW) 

=  n2bf  (1  + 

e  )  S(I,J) 
r 

Cf(IV)  = 

-i^nc  C(I ) 

C  (IW)  = 

-ipcf  jl(n2  ~ 

l)er  +  (n2  -  2)J  C( I ) 

where 

IU  = 

31-2 

JU  =  3J-2 

IV  = 

31-1 

JV  =  3J-1 

IW  = 

31 

JW  =  3J 

k  =  h2/12R2 

X  =  ttR/L 

b  =  Eh/R(l  -  v2)  t 

s 

:  =  Vs' 

S  =  W1 

-  v2)/hLR2 

e  =  e/R  b£  =  -EA^w  N  /R2L 

r  f  f  f  s 

Nf 

S(I,J)  = 

(2/N  )  I 

5  i=l 

sin  (Ittx^/L)  sin  (Jirx^/L) 

Nf 

C(I) 

(2/N  )  Z 

S  i=l 

sin  (Iwx./L)  sin  (I  irx./L) 

1  OC  1 

[A14] 


[A15] 


L^l] 


(Aid 


'oc  f 


[Aid 


[A20] 
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7.  Calculations  are  simplified  if  matrix  M  is  cast  in  the  form  of  equation  Q-0 

with  mUU  etc  defined  as  follows: 
o 


=  n2(l  +  k)  (1  -  v)/2  -  ♦n2bgSg 


uu  _  ,2 

m2  -  A 


uv 

mi 

= 

-nXU  +  v)/2  +  *  bsvAneg  -  nX/2j 

uw 

ml 

= 

vX  -  n2(l  -  v)kX/2  +  t|»X 

uw 

m3 

kX3 

w 

m 

o 

= 

2  „  -  sh 

n  -  *2n  bseQ 

vv 

m2 

= 

(1  -  v)(l  +  3k)X2/2  +  0sxX2bg 

vw 

m 

o 

= 

r  -  sh  sh] 

-n  ~  *n  |l  ~  2bseQ  +  bsVUx  J 

vw 

m2 

= 

-n(3  -  v)kX2/2 

ww 

m 

o 

= 

(1  +  k)  t  kn2(n2  -  2)  +  ^Jl  -  n  b 

ww 

m2 

= 

2  2  , 2- 

2n  kX  +  X  s 

S  O 

WW 

m4 

= 

kX4 

where : - 

iQ  =  aAsh(l  -  v2)/E 
9  9 

i  =  a  Sh(l  -  v2)/E 

x  x 


[A2f] 


m 


8.  With  reference  to  equations  [30]  and  {3l]  ,  the  c-vector  coefficients  approp 
riate  to  the  Kendrick  expressions  are  as  follows: 
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,  ,  -  sh 

-the  X  1  +  b  ve 
*  —  s  e 


oc 


s  e 
shl 


-tpe 


*eocnE  ’  Ve^ 

,<£  *  bsf<2  -  "2>;eSh  T  ^ 


-  l)u. 


/!] 


-  ,2- 

-I be  SAD 

OC  X  s 


JA23] 


w 


i^nc. 


~  o  2 

i he.  2(n  -  l)e  +  (n 

f  L.  r 


-«] 


CA24] 


V  wz  w  ,  , 

9.  Frame  coefficients  f  ,  f  and  f  are  given  by: 


rvw 


n2|f  +  ^2bjj 

-n|f{l  -  (n2  -  l)erl  +  #bf{2(l  +  er)  -  n2erj]  ^ 

f  jTl  +  e^)2  -  2n2e^  +  n2(n2  -  2)e^j  +  g(n  -  1)  +  ^bfn  (1  +  er) 
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